OU-HET 683/2010 

Entanglement Entropy of Two Black Holes and 
Entanglement Entropic Force 

Noburo ShibeQ 
Department of Physics, Graduate School of Science, 
Osaka University, Toyonaka, Osaka 560-0043, Japan 
(Dated: February 1, 2011) 

Abstract 

We study the entanglement entropy, Sc, of the massless free scalar field on the outside region C of 
two black holes A and B whose radii are R\ and R2 and how it depends on the distance, r(3> R\,R2) 
, between two black holes. If we can consider the entanglement entropy as thermodynamic entropy, 
' (— 1 ' ■ we can see the entropic force acting on the two black holes from the r dependence of Sc- We develop 

the computational method based on that of Bombelli et al to obtain the r dependence of Sc of 
scalar fields whose Lagrangian is quadratic with respect to the scalar fields. First we study Sc in 



d+ 1 dimensional Minkowski spacetime. In this case the state of the massless free scalar field is the 
Minkowski vacuum state and we replace two black holes by two imaginary spheres, and we take 



CO 
> 

o 

vo . 

, the trace over the degrees of freedom residing in the imaginary spheres. We obtain the leading 

CO ' 

term of Sc with respect to 1/r. The result is Sc = Sa + Sb + si-s G(Ri, R2), where Sa and Sb 
are the entanglement entropy on the inside region of A and B, and G{R\,R2) < 0. We do not 
calculate G(Rx,R<2) in detail, but we show how to calculate it. In the black hole case we use the 
method used in the Minkowski spacetime case with some modifications. We show that Sc can be 
expected to be the same form as that in the Minkowski spacetime case. But in the black hole case, 
Sa and Sb depend on r, so we do not fully obtain the r dependence of Sc- Finally we assume that 
the entanglement entropy can be regarded as thermodynamic entropy, and consider the entropic 
force acting on two black holes. We argue how to separate the entanglement entropic force from 
other force and how to cancel Sa and Sb whose r dependence are not obtained. Then we obtain 
the physical prediction which can be tested experimentally in principle. 

PACS numbers: 03.65.Ud, 04.70.Dy, 11.90. +t 
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I. INTRODUCTION 



Entanglement entropy in quantum field theory (QFT) was originally studied to explain 
the black hole entropy |2[ . Entanglement entropy is generally defined as the von Neumann 
entropy Sa = —TrpA In pa corresponding to the reduced density matrix pa of a subsystem 
A. When we consider quantum field theory in d + 1 dimensional spacetime M. x N, where 
R and N denote the time direction and the d dimensional space-like manifold respectively, 
we define the subsystem by a d dimensional domain A C iV at fixed time t = t . (So this 
is also called geometric entropy) Entanglement entropy naturally arises when we consider 
the black hole because we cannot obtain the information in the black hole. In fact, in the 
vacuum state the leading term of the entanglement entropy of A is proportional to the area 
of the boundary OA in many cases Q, 2I. This is similar to the black hole entropy, and 
extensive studies have been carried out 

In this paper we study the entanglement entropy, Sc, of the massless free scalar field 
on the outside region C of two black holes A and B whose radii are R\ and R2 and how 
it depends on the distance, r(^> R\,R2) , between two black holes. We consider the case 
that the state of the massless free scalar field is the vacuum state which depends how to 
choose the time coordinate. We choose the coordinate system which covers whole space 
time and does not have the coordinate singularity on the horizons. If we can consider the 
entanglement entropy as thermodynamic entropy, we can see the entropic force (we call this 
entanglement entropic force) acting on the two black holes from the r dependence of Sc- 

In Section [Til we obtain the general behavior of the entanglement entropy of two disjoint 
regions in translational invariant vacuum in general QFT. In Section UTTI we review the 
computational method of entanglement entropy in free scalar fie 



ds 111. There are some 



computation, method, of entanglement entropy |SHll| . See |l4 for reviews. That of 
Bombelli et al [1[ is most straightforward and powerful enough to obtain the r dependence 
of Sc in free scalar field theory. In Section [IV] we study Sc in d + 1 dimensional Minkowski 
spacetime. In this case the state of the massless free scalar field is the Minkowski vacuum 
state and we replace two black holes by two imaginary spheres, and we take the trace over the 
degrees of freedom residing in the imaginary spheres. We develop the method of Bombelli 
et al and obtain the leading term of Sc with respect to 1/r. The result in this section agrees 
with the general behavior in Section [TT] This method can be used for any scalar fields in 



curved space time whose Lagrangian is quadratic with respect to the scalar fields (i.e higher 
derivative terms can exist). In Section|V]we consider the black hole case. We use the method 
used in Section [IV] with some modifications. We show that Sc can be expected to be the 
same form as that in the Minkowski spacetime case. But in the black hole case Sa and Sb 
depend on r, so we do not fully obtain the r dependence of Sc- In Section ED we assume 
that the entanglement entropy can be regarded as thermodynamic entropy, and consider the 
entanglement entropic force. We argue how to separate the entanglement entropic force from 
other force and how to cancel Sa and Sb whose r dependence are not obtained. Then we 
obtain the physical prediction which can be tested experimentally in principle, and discuss 
the possibility to measure the entanglement entropic force. In Appendix B we obtain a 
formula for a finite series as a by-product of our calculation. 



II. GENERAL BEHAVIOR 



We consider entanglement entropy of two disjoint regions (A and B) in translational 
invariant vacuum in general QFT in d + 1 dimensional spacetime (d > 2). We will show 
that Sc reaches its maximum value when r — > oo. 

There are several useful properties which entanglement entropy enjoys generally. (See e.g. 
[ijj].) We summarize some of them for later use. 

1. If a composite system AB is in a pure state, then Sa = Sb- 

2. If pab = Pa® Pb, then S A b = S a + S B - 

3. For any subsystem A and B, the following inequalities hold: 

Sab < Sa + Sb, (1) 
Sab >\ Sa — Sb \ ■ (2) 

The first is the subadditivity inequality, and the second is the triangle inequality. 

Because of translational invariance, Sa and Sb are independent of their positions, 
so,^- = and ^f- = 0. And the total system is in a pure state, so we have Sc = Sab- 
Moreover, in the vacuum state, lim^oo pab = Pa® Pb because of the cluster decomposition 
principle [161 ] . So the property 2 suggests 

lim S c {r) = S A + S B - (3) 



r— >oo 
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We apply (GQ) and (j2J)to this system, then we obtain 

\S A -S B \<S c (r)<S A + S B . (4) 
Eqs. ([3]) and (j4j) show that 5c (as a function of r) reaches its maximum value when r — > oo. 



III. HOW TO COMPUTE ENTANGLEMENT ENTROPY 



In this section we review the computational method developed by Bombelli et al 



A. Entanglement entropy of a collection of coupled harmonic oscillators 

We model the scalar field on M. d as a collection of coupled oscillators on a lattice of space 
points, labeled by capital Latin indices, the displacement at each point giving the value of 
the scalar field there. In this case the Lagrangian can be given by 

L = \G MN q M q N ~ \v MN q M q N , (5) 

where q M gives the displacement of the Mth oscillator and q M its generalized velocity. The 
symmetric matrix Gmn is positive definite, and therefore invertible, i.e, there exists the 
inverse matrix G MN such that 

G MP G PN = 5 M N . (6) 

The matrix Vmn is also symmetric and positive definite. Next, consider the positive definite 
symmetric matrix Wmn defined by 

W MA G AB W BN = V MN . (7) 

The matrix W is the square root of V in the scalar product G. 

Now consider a region Q in M. d . The oscillators in this region will be specified by Greek 
letters, and those in the complement of Q will be specified by lowercase Latin letters. We 
will use the following notation 

, W ab W r\ I A B\ , r (\V ah W als \ ( D E\ 

Was = = W AB =[ = , (8) 

W ab W a J \B T C \W ab W a " \E T F 
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where W AB is the inverse matrix of Wab (W ab is not obtained by raising indices with G AB ). 
So we have 




A B\ D E\ / AD + BE T AE + BF 
B T C } \E T F] \ B T D + CE T B T E + CF 



(9) 



If we consider the information on the displacement of the oscillators inside Q as unavail- 
able, we can obtain a reduced density matrix p red for the outside Q, integrating out over R 
for each of the oscillators in the region Q, then we have 

Pre d (q a ,q' b ) = I n^VOf ,</WV), (10) 

J a 

where p is a density matrix of the total system. 

We can obtain the density matrix for the ground state by standard method, and it is a 
Gaussian density matrix. Then, p re d is obtained by a Gaussian integral , and it is also a 
Gaussian density matrix. The entanglement entropy, S = —trp re d\n p re d , is given by Q 

S = £/(A n ), (11) 

n 

/(A) = ln^A 1 ' 2 ) + (1 + A) 1/2 ln[(l + A" 1 ) 1 / 2 + A" 1 / 2 ], (12) 

where A n are the eigenvalues of the matrix 

A a b = -W aa W ab = —(EB T ) a b . (13) 

It can be shown that all of A n are nonnegative as follows. From OH]) we have 

AA = —AEB T = BFB T . (14) 

It is easy to show that A, C, D and F are positive definite matrices when W and W~ x are 
positive definite matrices. Then AA is a positive semi definite matrix as can be seen from 
(I14p . So all eigenvalues of A are nonnegative. Finally, we can obtain the entanglement 
entropy by solving the eigenvalue problem of A. 



B. The continuum limit 



Next, we apply the above formalism to a massless free scalar field in (d+1) dimensional 
Minkowski spacetime. We take the continuum limit in the above formalism. In this case the 
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Lagrangian is given by 

L = J d d x\[tf - (V0) 2 ]. (15) 
Then the potential term becomes 

\v AB q A q B -> J d d x l -[{Vct>f]. (16) 
The matrices V, W and W^ -1 are given in the momentum representation by, 

W(x,y) = J 0- d (k 2 ) 1/2 e ik <^ (18) 

W~\x,y) = J ^ d (kr i/2 e^ (19) 

From (1131) . the matrix A is obtained as a sum over the oscillators in the region Q, 

A(x,y) = - I d d zW- 1 {x,z)W{z,y). (20) 
Jq 

We now have to solve the eigenvalue equation 

/ d*yA(x,y)f(y) = \f(x), (21) 

where Q c is the complementary set of Q, and then we use the eigenvalues in the expression 
for the entropy ( ITT]) . 

IV. ENTANGLEMENT ENTROPY OF TWO DISJOINT REGIONS IN A d + 1 
DIMENSIONAL MASSLESS FREE SCALAR FIELD 

We consider two spheres A and B whose radii are R\ and R2, and define the outside 
region as C. (See Fig CD) We derive the r(^> i?i,i? 2 ) dependence of Sc(r, Ri, R2) by using 
the formalism of the preceding section. (In the later analysis we do not use the shapes of 
A and B, so all analysis in this section holds for A and B which have arbitrary shapes. In 
this case R\ and R2 are the characteristic sizes of A and B.) 

We consider Sab^, Ri, R2) because Sc = Sab in a pure state and the r dependence of 
Sab is clearer than that of Sc in the calculation. In this case the region Q is C. 

We obtain the r dependence of Sab by following three steps: 
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(1) We obtain the r dependence of A by using the ||x — y\\ dependences of W(x, y) and 
W /_1 (a;, y). We decompose A into the non-perturbative part and the perturbative part as 
A = A + SA, where A = lim^^ A. 

(2) We obtain A m (r) which are the eigenvalues of A by perturbation theory. This is almost 
similar to the time-independent perturbation theory in quantum mechanics in the presence 
of degeneracy. We can regard A as Hamiltonian. Note that A is not a symmetric matrix. 
So we must slightly modify the perturbation theory in quantum mechanics. 

(3) In Step (2), we had A m (r) as A m (r) = A^ + 5X m (r), where A^ are the eigenvalues of 
A . We substitute these A m (r) into (TT21) . then we obtain Sab^, Ri, Rz)- 

First we examine the — y\\ dependences of W(x^y) and W~ 1 (x,y). Generally en- 
tanglement entropy has UV divergence as discussed in l). So we use a momentum cutoff 
l~ x in integrals (jT71)-(|T9l). though these integrals are well defined as Fourier transforms of 
distributions. (The other regularization methods are discussed in When d > 2 and 
l/\\x — y\\ —7-0 , W(x, y) and W~ l {x, y) are 



W(x,y)= J 4 ' , W-\x,y)= Bd A d ,B d eR (22) 

\\x y\\ \\x y\\ 



where A d and B d are nonzero dimensionless constants (see Appendix A). We cannot ob- 
tain ( )22|) by only using a dimensional analysis because l/\\x — y\\ is dimensionless. Indeed 
V(x,y) — > when l/\\x — y\\ — > 0, i.e. V(x,y) is zero when ||x — y\\ is finite. On the other 
hand W and W~ x have nonzero value for ||x — y\ > because they are kernels of integral 
operators of nonlocal interaction (i.e Fourier transformations of (y/k?)^ 1 ) . In Appendix A 
we explicitly show that W and W~ x have nonzero value for — y\\ > and Eq. (1221) holds. 

Next, we obtain the r dependence of A by using ff22l . We represent the matrix elements 
of iy(iy _1 ) diagrammatically in Fig El Instead of solid lines we use dotted lines for W~ l . 
The lines denote the matrix elements W(x, y)(or W~ 1 (x, y)) in ( 1221) . An initial point and an 
end point of an arrow denote a row and a column respectively. We can obtain products of 
matrices by connecting arrows and integrating joint points on regions where the joint points 
exist. We label coordinates in A, B and C as x a ,Xb and x c . Then, from Fig [3] we obtain 





c 






) 
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FIG. 1. Two spheres A and B, and the outside region C. 
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FIG. 2. The matrix elements of W. The lines denote the matrix elements W(x,y) in (|22h . An 
initial point and an end point of an arrow denote a row and a column respectively. We can obtain 
products of matrices by connecting arrows and integrating joint points on regions where the joint 
points exist. Instead of solid lines we use dotted lines for W -1 . 

A = —EB T as 



A 



W{z c ,y a ) W{z c ,y b ) 



(23) 



A(x a ,y a ) A(x a ,y b )\ (W 1 (x a ,z cj 
A(x b ,y a ) A(x b ,y b )J \W~ 1 (x b ,z c 

f c d d z c W~ 1 {x a , z c )W(z c , y a ) f c d d z c W~ 1 (x a , z c )W(z c , y b ) 
f c d d z c W-\x b , z c )W(z Cl y a ) J c d d z c W-\x b , z c )W{z C) y b ) 

To make the r dependence of the non-diagonal elements of A clear, we use the following 
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A = 
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FIG. 3. The diagrammatic calculation of A in ([23 
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FIG. 4. The diagrammatic representations of A in (|25p and the identity in 
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A = A D + SA t + 5A 2 



FIG. 5. The diagrammatic representations of <5Ai, 5A2 and Ap. 
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FIG. 6. The diagrammatic representations of A , A and A°5Ai. 
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FIG. 7. The diagrammatic calculation of 5 Ad in (|34p . 



identity, 

/ d d zW-\x a , z)W(z, y b ) = 8{x a - y b ) = 0. (24) 

J A+B+C 

We represent this identity diagrammatically in Fig HI From ( 123|) and ( 1241) we obtain (see 
Figi]) 

A(x a ,y b ) = / d d z a W-\x a ,z a )W(z a ,y b ) + / d d z b W~ 1 {x a , z b )W{z b ,y b ) 

Ja Jb (25) 

A(x b ,y a )= / d^a^- 1 ^,^)^^,^) + / d d z b W- 1 (x b ,z b )W(z b ,y a ). 
J a Jb 

Note that from ( )22|) PF(x,?/) and VF _1 (:r, y) have the different — y|| dependence. So, from 
and f l23|) we decompose A as 



A = A D + 5A 1 + 5A 2 
10 



(26) 



where we define (see Fig [5]) 

A D =i AM ° ], (27) 
A(x h ,y b )^ 

M , j B d d z b W-\x a ,z b )W(z b ,y b ) 

j A d d z a W-\x b ,z a )W{z a) y a ) 

M , f A d d z a W-\x a ,z a )W(z a ,y b ) (2Q) 

J B d d z b W~\x bl z b )W(z bl y a ) 

We approximate W(x a ,y b ) ~ -3+r an d W 7-1 ^, 2/fc) ~ ^hr because r 3> Rx,R 2 . Then we 
have 

\S A d d z«W(z a ,y«) 



r" 



5A 2 « 4r7 . 31 

Next we consider the non-perturbative part A = lim^oo A. From ( 1301) and ( 13T1) we can 
see that 5Ai and 5A2 become when r — » 00. Note that the integral region of the integral 
in A(x a , y a ) (A(xb, 2/5)) become A c = M. d — A (B c = M. d — B) when r — > 00, then we obtain 
(see Fig [6]) 

A o = _ (f Ac d d zW-\x a ,z)W(z } y a ) 

/^c^w-V^W^) 



From ( 1321) we rewrite ( |26l) as follows, 

A = A + 5Ai + 5A 2 + 5A D , (33) 

where we define (see Fig [TJ) 

/ f R d d z b W- 1 (x a ,z b )W(z b ,y a ) 
5A D = A D - A = \ Jb y ! V y J . (34) 

\ J A d d z a W-\x b ,z a )W(z a ,y b ) 

We use the same approximation as we used in ( 130]) and ( I3"T1) . then we obtain 



f A c^ a 
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When we perform the perturbative calculation to obtain A m (r) which is the eigenvalues of 
A, from fl30l) . (l3Tj) . f )33|) and fl35|) we can neglect because it is higher order than 5Ai 
and 5A 2 with respect to 1/r. And we can neglect 5A 2 because its nonzero matrix elements 
are in the same position as 5Ai and 5A 2 is higher order than 5Ai with respect to 1/r. 

Because A is not a symmetric matrix, in the later perturbative calculation we need A°5Ai 
where A is defined as (see Fig [6]) 



A = lim A 



W{x a ,y a ) 







W(x b ,y b ) 



(36) 



From (J2HD, flU and (J36J we obtain (see Fig[BJ 





J A (i d z a ^^^(aJa,^)^ 1 (z a ,z b )W(z b ,y b ) 



jB ddz b J A d d z a W(x b ,z b )W 1 (z b ,z a )W(z a ,y a ) 



B d 



J A d d z a W(x a ,z a ) J B d d z b W(z b ,y b ) 



J B d d z b W(x b ,z b ) f A d d z a W{z 



(37) 



We have finished the first step. 

Next, we calculate A m (r) by perturbation theory. This is almost similar to the time- 
independent perturbation theory with degeneracy of quantum mechanics. The only differ- 
ence is that A is not a symmetric matrix and AA is a symmetric matrix. 

We can approximate A rj A + SAi and regard 5Ai as the perturbative part. Then, from 
fl30|) we expand A m (r) with respect to l/r d_1 . We expand A m around A^ = A m (r = oo), 



A m — A^ + 5A^ + 5X 2 m) 



(38) 



where 5A^ and <5A^ are the first and the second order perturbations. 
Next we substitute ([38} into (PT2D, 



S , Afl(r,i2 1 ,ife) = 53/(A w ») 

m 

= 5 il (i2i) + 5 , B(^2) + Z) 



<5A r 



dX r 



Am — A™ 



+ 7^A m ) 2 



dAJ 



Am — A? r , 



(39) 



where 5X m = 5A^ + 5A m . We will show that the first order perturbation in (139]) (i.e 
' ) is zero, so we must calculate the second order perturbations. 

Am, — A™ 



J2 m S^m 
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We label the A^'s as A^ > A° when m > n. And we define the eigenvectors of A , 

-0 _ (fmla( x a)\ [ \ _ ^ 

Jmla ~~ I I iJn2l3 — I \ I ' ■> mla ~ A mJ mlai lv J m2/3 ~ A mJ m2fii 

where a — 1, • • • , M m and /3 = 1, ■ ■ • , N n are the labels of the degeneracy. And we normalize 
fLa (i = 1,2) as follow, 

f%«A &f> = tmJiAe- (4i) 

This normalization is always possible because A is a positive definite symmetric matrix. 
For general Ri and R2, A 01 and A 02 have different eigenvalues, so there are two groups of 
A^; one is the group of the common eigenvalues of A 01 and A 02 , the other is not. We will 
see that 5A^ of the latter group are zero. We expand / m7 which is the eigenvector of A in 
the following way, 

fmy — ^ ] a yaf m la ^ ] ^pfmSP fmy fmy = £7717 fmy frrvy (^2) 
a 

where and / 2 are the first and the second order perturbations. Note that when A^ 
is an eigenvalue of A 01 (A 02 ) and is not an eigenvalue of A 02 (A 01 ), then the coefficients 6 7j g 
(a 7Q ) are zero ; because the zeroth order eigenvectors f^ 2/3 (f m i a ) do not exist. So either 
the coefficients a ia or b 1 p are zero when A^ is not a common eigenvalue of A 01 and A 02 . We 
substitute (H2l) into the eigenvalue equation (we approximate A rs A + SA{) , then we have 

(A + 5A 1 )f nn = {X° m + SX 1 ^ + 5\ 2 mi )f mi . (43) 

We obtain equations of the first and the second order perturbation. 

AO fl I n t0 _ \Q el _|_ r\l tQ 

J my * lsm7 mJ my ' my^my \ J 

J my ' " lim7 A mJmy ' ^ A myi my ' ^ A my^>my v*"/ 

We multiply fj44l) by ,A° from the left . The first term of the left hand side of (jHJ) cancel 
the first term of the right hand side of (jHJ) because A A is a symmetric matrix, then we 
obtain 

E + E ^Kf ym/3 = SX^ia^ 1 + 6 77 , <P 2 ), (46) 

a $ 
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where 



yij — fOT ^Oca fO 



(47) 



From (J37D we obtain V^ an/3 - v manp 



V 



22 



and 



(^fmla( X a) ®j 



A°5A 1 (x a ,y a ) A°5Ai(x a ,y b ) 
A°5Ai(x b ,y a ) A°5Ai(x 6 ,y 6 ) ) \f^(y b ) 



= -zk I ddx a I d d z a W(x a ,z a )f^ la (x a ) [ d d y b [ d d z b W {y b) z b ) fl w {y b ) = -^-C man/3 

r J A J A JB JB r 

(48) 

and V^ anj3 = V^ ma . We define an M m x N n matrix C mn as [C mn ) afi = C man p and write 
f|46|) as follows, 



r d-i 1 „ T 




(49) 



b 7 / \ b 7 

where (a 7 ) a = a 7Q , and (h y )p = 6 7/ g. From ( 1491) . if A^ is not a common eigenvalue of A 01 and 
A 02 , 5A* is zero; because either a ia or 6 7j g are zero when A^ is not a common eigenvalue of 
A 01 and A 02 . We first consider the case that M m > N m . In this case we obtain the following 



eigenvalue equation 



det 



xl 



Af m X JV/ m 



-C T rl 



N m xN m 



det(a;l 



M m xM m 



det(xl 



N m xN m 



(50) 



„M m -N m 



X 



det(x l NmXNm - C mm C mm ) - 0. 

We define the eigenvalues of C^ m C mm as c ma (a = 1, ■ ■ ■ M m ). C^ im C mm is a positive 
semidefinite matrix because C mm is a real matrix, so c ma > 0. Then we obtain 5 A* from 



91) and (1501) ■ 



5Xi 







rr»7 



5, 



(51) 



±—[Zi^fc^ (o£=l,---M m ) 

When M m < N m , we can obtain 5A^ 7 in the same way. We define the eigenvalues of 
C mm C^ m as d ma (> 0) (a = 1, ■ • • iV m ). Then we obtain 



5Xi 







j«7 



5, 



(52) 



fa=l,---JV„ 
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Then £ *A 



x df 



m )7 m 7 d . 



Am — A m 



0, because we have ^ 7 £A^ 7 = from (T5TT) and 

Next we consider £A^ 7 . We skip the detailed calculation because it is also almost similar 
to the time-independent perturbation theory with degeneracy of quantum mechanics. Then 



we can write £A^ 7 as follows 



X\2 



E 



\0 _ \0 

n{^m),i,p m n 



(/nT^°5AiC 7 )(e^°5Ai/^) 



m n 



(53) 



where 



— E/ fnififnlpA ■ 



(54) 



n is a projection operator on the eigenspace of A°. To obtain <5A^ 7 we must obtain £° by 
solving the eigenvalue problem, but it is not necessary for our purpose because we want to 

. From (153]) we obtain 



> d f 



know only £ m , 7 5A m 7 rfA 



A m =A« 



m, 7 



2 

m7 ciA 



E 



m,n(mj^n) 



A^-A° 



c/A r 



E 



A - A 

m,n(rn>n) m n 



Tr((f) m SA 1 4) n SA 1 



1 df 








\ d\ m 


^?7] — -^m 


d\ n 





In the second line we have used 



^m 7 ^m 7 



(55) 



(56) 



and in the third line we have used cyclic property of trace. Next we examine the sign of 
( )55l) . Its trace term is positive because 



Tr^SA^JA,) = E (C^SAJ^ifJ^A^AX, 

= y v ij R v- 



' m/3na / j \ nam/3 



2 > 0. 



(57) 
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And from (1121) we obtain 
df 1 



d\ 2V1TX 

ll_ 1_ 

d\ 2 ' 



111 



1 1 

1 + T + ^f 



> for A > 0, 



4V1 + A 



1 + A 



In 



1 1 

l + T + ^= 



+ 



AvO + A 



(58) 

< for A > 0. (59) 



From ( ]5"7|) . (159]) and A^ > A° (m > n) , (155]) is negative. And from ( 159]) we obtain 

<Pf 



m,7 



™ 7 dA 2 



< 0. 



(60) 



Am — A^, 



Finally, from (]3"9"]) ; AST]), (152D, (153]) and (1ST]) we obtain 

df 



S AB (r, R u R 2 ) - - S B (i? 2 ) = ^ 



m,7 



X\2 

m7 dA r 



A m — A^. 



A m =A0 



5, 



„<i-l 



£(o, 



^0 _ yi £_^\~»<<"< 

m,n(m>n) m n a, (3 



I df 




d/ 




\ dXrn 




d\ n 


An=A0/ 



(61) 



+ £ 



Cm 1 r 



d 2 f 



d\ 2 m 



m"a 



d 2 f 



Am=A°, m",o m 



1 



-2d-2 



G{R U R 2 ) < 



where denotes the summation taken over the common eigenvalues of A 01 and A 02 , 

whose degeneracy is M m > N m , and Ylim" denotes the summation taken over the common 
eigenvalues of A 01 and A 02 , whose degeneracy is M m < N m . 

We have obtained the r dependence of Sc(r, R±, R 2 ) = Sab^, Ri, R2) m ( 161]) . then we 
next consider G(R\,R 2 ). To calculate G(R\,R 2 ) we need to know A^ and f^ ia which we 
do not examine in this paper. But from C man p{R\ = 0,R 2 ) = C man p(Ri, R 2 = 0) = we 
obtain a trivial property of G(Ri, R 2 ), 



G{Rr = 0, R 2 ) = G{R±, R 2 = 0) = 0. 



(62) 



And G(Ri,R 2 ) depends on the cutoff length I because A^ , f^ ia and C man/ 3 depend on I. 
(A^ are dimensionless, so they depend on R\/l or R 2 /l. And in (115]) f A d d z a W(x a ,z a ) and 
f B d d ZbW(yb, Zb) depend on / because W(x,y) depend on I for x y, so C mQ , n/ 3 depends on 
/. ) Probably G(R±, R 2 ) diverges when I — > 0, as Sa(Ri) and Sb(R2) have divergence 
[l, 2|. And G(Ri,R 2 ) most likely diverges more weakly than S^C^i) an d ^(-^2)- Then, by 
dimensional analysis, when Ri = R 2 = R we can assume 



G(jRl = i2, # 2 = fl ) = ^ 2 ^ 2 ( ^ m An ^ 



d-l>m>0,n>0,g<0 (63) 
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FIG. 8. The diagrammatic representations of Ad, A and 5 Ad- 

where g is a dimensionless constant. 

Finally we consider the condition under which the approximations are good. When 



r ^> i?i,i?2, SA « 5Ai is a good approximation. When 



B d 

r.d— 1 



man/3 



< |Ai - A°|, the 



perturbation theory is a good approximation. The latter condition might have I dependence, 
so we might need the condition R/r <C (l/R) a , where a > 0. 

V. ENTANGLEMENT ENTROPY OF TWO BLACK HOLES IN A d + 1 DIMEN- 
SIONAL MASSLESS FREE SCALAR FIELD 

In this section we consider the entanglement entropy of the massless free scalar field on 
the outside region C of two black holes A and B whose radii are Ri and i?2- The action of 
the massless free scalar field is given by 



S 



(64) 



First we specify the vacuum state of the scalar field. The vacuum state is specified by spec- 
ifying the time coordinate t. We use the coordinate system which have following properties: 
this coordinate system covers the inside and the outside regions of two black holes and does 
not have the coordinate singularity on the horizons and becomes the orthogonal coordinate 
system of Minkowski spacetime in the region far from the two black holes. To construct 
this coordinate system, we use the coordinates which is similar to the Kruskal coordinates 
in the inside regions and the neighborhood of black holes, and similar to the Schwarzschild 
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coordinates in the other region. In this coordinate system g a is positive everywhere, then 
from (jMj) Gmn and Vmn in © are positive definite. So we can use the formalism in the 
Section IHIl 

We can use the method of the last section with some modifications. In the black hole 
case W(x, y) and W~ 1 (x, y) depend on r, so we write them as W(x, y;r) and W~ l (x,y;r). 
Exactly in the same way as in Minkowski spacetime, Eqs. ( 126|) - (|29|) hold because (124)) holds. 
On the other hand Ao(= lim^oo A) changes because W(x, y; r) and W~ 1 (x, y; r) depend on 
r. We define WA(x,y) and (x,y)(WB{x,y) and W B (x,y)) as W(x, y) and 
in the case that the only one black hole A(B) exists. Then we have 

A0 = _ff Ac d d zWX 1 (x a ,z)W A (z,y a ) \ 

\ J Bc d d zW B 1 (x b ,z)W B (z,y b )) ' 

It is difficult to evaluate the r dependence of 5Ae> = A^ — A because it is difficult to evaluate 
W(x,y;r) — WA(B)(x,y) and W~ l (x, y;r) — W^ B ^(x,y). So, in the black hole case we do 
not consider Aq as the non-perturbative part. Instead we define Ad(t) and A(r) as (see Fig 



A D {r) 



A A {x a ,y a ;r) 

A B (x b ,y b ;r) 

-J Ac d d zW- 1 {x a ,z;r)W{z,y a ,r) 

- J Bc d d zW- 1 (x b ,z;r)W(z,y b ;r) 



(66) 



(W(x a ,y a :r) 

A{r) = 1 ' y ' ; | , (67) 

y W(x b ,y b ;r)^ 

and we consider Ap. as the non-perturbative part. Note that A^ and A# are the matrices 
A corresponding to Sa{t, Ri, R2) and Sb{t, Ri, i^)- So we will obtain Sab as the following 
form, Sab{t, Ri, R2) = Sa(j, R\, R2) + Sb(t, Ri, R2) + o~Sab{j~, Ri, -^2)- We calculate the 
leading term of SSab^, Ri, R2) with respect to 1/r. 

We define SA D = A D - A D (see Fig ED, then we have A = A D + 8A t + 5A 2 + SA D . To 
evaluate 5Ai, 5A 2 and 5A D , we evaluate W(x a , y b ; r) and W /_1 (x a , y b ; r). When r 3> Ri, R 2 , 



by dimensional analysis we obtain W(x a ,y b ;r) - 3 ^ r L 1 (i? 1 /r, R 2 /r) and 1 (x a ,y b ; 



r. 



, . 3 - 



jL2{R\/r,R2/r), where Li and L2 are dimensionless functions of -Ri/r and R2IT. The 



space time becomes Minkowski space time when Ri — y and -R2 — ► 0, so in this limit 
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A8A 1 



yd— 1 



probably we have L\ — > 1 and L2 — >■ 1. This limit is equivalent to r — > 00, so we have 
lim^ooLi = lim^ooZ^ = 1- Then we obtain 8A 1 = OiA/r 11 ' 1 ), 5A 2 = 0{l/r d+1 ) and 
5 = 0(1 /r 2d ) as well as the Minkowski spacetime case. We can neglect 5A 2 and 5Ad 
for the same reason as in the Minkowski spacetime case (see below Eq. (|35p ). So we can 
approximate A »s Ap. + SA\. Then we change the perturbative calculation in the last section 
as follow 

A°^A D (r) A°^A(r) X° m ->■ X° m {r) f mia ^ /° ia (r) (68) 

where A^(r) and fmia( r ) are the eigenvalues and the eigenvectors of A D (r). The perturbative 
calculation is the same as that in the last section. In this case A(r), A^(r) and fmia( r ) depend 
on r, but we can remove their r dependence as follow. Because we want to calculate the 
leading term of S AB (r, Ri, R2) — S A (r, R 1 , R 2 ) — S B (r, R 1 , R 2 ) with respect to 1/r, we can 
approximate 

B d L 2 (f,f) I J A d d z a W(x a ,z a ;r) J B d d z b W(z b ,y b ;r) 

J B d d z b W(x b , z b ; r) J A d d z a W(z a , y a \ r) 

B d I f A d d z a W A (x a ,z a ) J B d d z b W B (z b ,y b ) 

rd ' 1 \f B d d z b W B (x b ,z b ) f A d d z a W A ( y z a ,y a ) 

(69) 

In the second line we have approximated L 2 (^-, m 1, W(x a , z a ;r) rs W A (x a ,z a ) and 
W(z b , y b ; r) « W B (z b ,y b ). And we can approximate A^(r) w A^(r = 00) = A^ and 

/U r ) ~ 7mia( r = 00 ) = /la- Note that A m and /mi« are the eigenvalues and the 
eigenvectors of A , i.e. (A is in ( |65l) ) 

„0 _ I fmla( X a)\ rO _ ( ® | aOj-0 _ \0 fO A f — \° f° f7n"\ 

imla ~~ I I iJn.2/3 ~ I \ I ' -/mla ~~ A mJ ml«> iv Jm2/3 ~~ A mJ m2(3i V U J 

\ / \fn2p( X b)J 

where a = 1, ■ • • , M m and = 1, ■ • • , N n are the labels of the degeneracy. 
Finally we obtain 

S AB (r, R ± , R 2 ) = S A (r, R u R 2 ) + S B (r, R u R 2 ) + ^G(R U R 2 ) (71) 

where G(R\, R 2 ) is the same function as that in (IBTj) . Note that in this case from ( 169]) C man p 
in G(Ri, R 2 ) is 

C manP = / d d x a / cf*z W A (z , z a )f^ la {x a ) / / ^ 6 Wb(2/6,^)/°2^(2/6)- (72) 
J A Ja J_b 
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As in the Minkowski spacetime case, we obtain G(Ri = 0, R2) = G(R\,R2 = 0) = from 
CmanpiRi = 0, -R2) = C mC mji{R\-, Ri — 0) = 0, and G(Ri, R2) probably diverges when I — > 0, 
where / is the cutoff length. The l/l dependence of G(Ri,R 2 ) is most likely the same as 
that in the Minkowski spacetime, then we obtain 

G(Ri = R ) R 2 = R)= g BH R 2d - 2 (^) fin ( *f)) d-l>m>0,n> 0,g BH < (73) 



j j \ v 1 , , 

where qbh is a dimensionless constant, and m and n are the same numbers as those in the 
Minkowski spacetime. 

VI. ENTANGLEMENT ENTROPIC FORCE AND THE PHYSICAL PREDIC- 
TION 

We assume that we can consider the entanglement entropy of two black holes as thermo- 
dynamic entropy. If this assumption is correct, the entropic force acts on two black holes. 
We consider the force of the scalar field which acts on two black holes. We consider two 
black holes which have same radius R\ = R2 = R, then we can consider the temperature T 
to be the Hawking temperature. We define the energy and the free energy of the field on 
the region C as Ec(r, R) and Fc(r, R), 

F c (r, R) = E c (r, R) - TS c (r, R) = E c (r, R)-T \2S A {r, R) + ^G(R)^j . (74) 

where G(R) = G(Ri = R,R 2 = R) and we have used (1711) . We define the force of the field 
on the region C which acts on one black hole in the direction of increasing r as Xc- We 
obtain Xc by partially differentiating Fq with R fixed, 



(75) 



In (175|) the second term is the entropic force. 

We cannot see the effect of the entropic force only from (1751) because we do not know 
Sa(t,R). To see the effect of the entropic force we consider three situations . (See Fig[H]) 
(1) There are two black holes which have the same radius R and the distance between them 
is r.(This is the situation we have considered.) (2) There are one black hole whose radius 
is R and one solid ball whose radius is Rq ~ -R(-Ro > R)> an d the distance between them is 
r. This ball has mass M which is the same as that of a black hole whose radius is R. And 
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FIG. 9. Three situations to see the effect of the entropic force. (1) There are two black holes. (2) 
There are one black hole and one solid ball. (3) There are two solid balls. We define the force of 
the field which acts on one black hole or on one ball in the direction of increasing r as , X^ 
and 



the scalar field does not exist in this ball. The boundary condition on the scalar field on 
the surface of this ball is not so important in the later calculation that we do not specify 
the boundary condition in detail. We only require that the scalar field on the outside region 
of this ball is not so different from that in the situation (1). (3) There are two solid balls 
which have the same radius Rq and the distance between them is r. These balls have the 
same properties as those in the situation (2). 

We define the force of the field which acts on one black hole or on one ball in the direction 
of increasing r as Xq\ X(9 and X^. We illustrate in Fig [9] the directions of force and the 
names of the regions. 

In the situation (2) the state of the field is \0) A+C2 , where [0)^^ is the vacuum state 
on A + C 2 . Because \0)^ic 2 

is a pure state, then = Sf. We define A^Aj?) as A 
corresponding to (S^)- Because the scalar field does not exist in the ball, then we 
obtain 



A? -A? 



^d d z b W-\x a ,z b )W(z b ,y a ) = O . (76) 
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Then we can approximate = S^\r,R) + 0(-^) ~ S A \r, R). Then we obtain 
X g\r,R)- - 9E ^ R) | T dS A^R) ^ ggj^ +T ggj^ _ (?7) 



In the situation (3) the state of the field on the region C 3 is a pure state, so = 0. 
Then we obtain 

(3) dF® dE^(r,R) 
Xg(r,i2) = ~^ = (78) 

From ([7SD dZZD and (EHD we obtain 



Xg> - 2X% + X® « -^[4? - 24? + 4i] - (2d - 2)T-^G(R). (79) 

o P (2) , P (3) 



— 2E , £. 2 + is Casimir energy. 
We have not considered the force of gravity. But we can include them in (IT9l easily. We 
define total force acting on one black hole or on one ball in the direction of increasing r as 
Ja\ J- a and F^n- Then we obtain 

jf-2J$ ) +J^=X$-2Xg>-X$> « -| r [4 1 i ) -24)+4)]-(2d-2)T-i rT G(/?). 

(80) 

The force of gravity is canceled in (IHUl) . The first and the second terms in the right hand 
side are the Casimir force and the effect of entropic force, respectively. 



Finally we consider the case d = 3. In this case the Hawking temperature is T 
From (J73]) and (ED]) we obtain 



1 



8nG N M 



o F (2) , F (3) _ 5 rpW _9P (2) 4- F (3) l - 9bhR 3 (R\ m f, (R\\ 
T A - 27 A + ^ ~ -— [E Ci 2E C2 + | __ X \T ) ) 



(81) 



2 > m > 0,n > 0,^bh < 0. 



We roughly estimate the Casimir force by analogy with that of electromagnetic field between 
two dielectric spheres with center-to-center distance r in Minkowski spacetime. The Casimir 
force between the two sphere was calculated in [15j, and it is 0(l/r 8 ). So, in our case we can 
probably neglect —§ ? [E§ ) 1 -2E i £ + E^} in fl8TJ). The left hand side of QgQ can be measured 
experimentally, so ( 1ST]) is the physical prediction. From ( 1ST]) the effect of the entropic force 
becomes significant when R is large. We can probably use heavy stars as the balls in the 
situation (2) and (3). So we can possibly confirm the effect of the entropic force by the 
cosmic observation (e.g. binary black holes and binary neutron stars). 
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We estimate the magnitude of the effect of the entropic force. We set the cutoff length I 
to the Planck length lp = (Gat/j/c 3 ) 1 / 2 , then the ratio of the effect of the entropic force to 
the force of gravity is 

Feef _ ^9bh (lp) 2 R ( R\ m ( ' f R 



Mr (82) 



T g n r 3 \l P J \ \l P 



where 



7T r 5 VW V. VpJ J r 2 4G N r< 

If m = 2, Zp is canceled in f l82|) . In this case the effect of the entropic force is comparable to 
the force of gravity, so we can possibly observe the effect of the entropic force. 



VII. CONCLUSION AND DISCUSSION 

In Section[TTlwe showed that the entanglement entropy (So = Sab) of two disjoint regions 
in translational invariant vacuum in general QFT reaches its maximum value when r — > oo. 
And we obtained the inequality (j3J). In Section [IV] we developed the method to obtain the r 
dependence of Sc and obtained the r dependence of Sc (1611) in the free massless scalar field 
in (d + 1) dimensional Minkowski spacetime. We can use this method in curved space time 
and for scalar field theory whose Lagrangian is quadratic. To know only the r dependence we 
need only the — y|| dependence of W(x,y) and iy _1 (x,y) when \\x — y\\ is large. To know 
the Ri and i?2 dependence we must solve the zeroth order eigenvalue equation and obtain 
}P m and fmia- ft i s difficult to solve the zeroth order eigenvalue equation analytically so we 
will need to perform numerical calculation. But we assumed the R\ and R2 dependence ( 1631) 
by using dimensional analysis and the cutoff dependence of Sa and Sp- 
in Section|V]we showed that Sc can be expected to be the form (TfB in the black hole case. 
In this case the only assumption we made is the r dependence of W(x a , yi) and iy _1 (x a , ?/&). 
We did not explicitly calculate W(x a , y^) and W~ l (x a , y^), but assumed the r dependence 
of W(x a ,yb) and W /_1 (x a ,?/b) by dimensional analysis. 

In Section IVII we assumed that we can consider the entanglement entropy of two black 
holes as thermodynamic entropy, and investigated its entropic force. We considered three 
situations (1), (2) and (3) and obtain the relationship (I8ip between the force acting on one 
black hole or on one ball and the sum of the Casimir force and the effect of the entanglement 
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entropic force. Because we can probably neglect the Casimir force, we can confirm (IHTj) 
experimentally in principle. And we can possibly confirm the effect of the entropic force by 
the cosmic observation because it is significant for large black holes. 

Next we discuss the entanglement entropic force in different systems. In the black hole 
case, black holes act as "walls" which hide inside regions but hold the entanglement between 
inside and outside regions. So if there are walls of this type, the entanglement entropic 
force will exist between regions surrounded by these walls. Then we will be able to confirm 
the entanglement entropic force by experiments in a laboratory if we make this wall. And 
if entanglement entropy depends on some external parameter, entanglement entropic force 
probably appears also in quantum mechanical (i.e. not quantum field theoretical) systems. 

Finally we mention our assumption that we can consider the entanglement entropy of two 
black holes as thermodynamic entropy. Entanglement entropy has property which is different 
from that of thermodynamic entropy. For example entanglement entropy is not a extensive 
variable in general. So we must reconsider statistical mechanics from a fundamental level to 
judge whether our assumption is correct or not. We can also use (1ST]) to judge the correctness 
of our assumption by experiments. 
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Appendix A: The calculation of W and W~ x 

In this appendix we calculate W(x,y) and W~ l {x,y) ( ffTBj) and f[T§j) ) explicitly. We 
regularize them by including convergence factor e~ 1 ^ in them, where I is the cutoff length. 
We define W a as 




(Al) 



Then we have Wq = W and W2 = W x . First we consider the case d > 3. 
(1) d> 3 



24 



-iz 



C 





1 f 


-1 






r-r ► 

J 

Cr 



{a) 



(b) 



FIG. 10. The contours of the integrals, (a) d = 2m + 2 (m > 1). (b) d = 2m + 1 (m > 1). 

We perform the integrals of angular coordinates which do not enter the inner product, 

/ 



W a (x,y) 



(2tt)< 



d-2 

n 

m-2 



r 



V 



r 



d — \ 
/ 



cZ — m + 1 



/l ^— 3 

dt[l - t 2 ]~e ifcr W" Q (A2) 
i 



where r = \\x — y\\ and we change the variable as t = cos8. Next we perform the k integral 



oo pi d _ 3 rl 

dk / dt[l - t 2 }~e ikrt e~ lk k d - a = 



o 



,.,-,,-r, j_ ■ (t + iz) d ~ a+1 



= H) d - a -\d-a)\ ?r ^ I dt\l-t 2 ]~ 
where z = Z/r. We define 



dt\l-t 2 }~ . 
i J \itdr 

d-3 



d—ct 



ikrt —Ik 



dke lkrt e 



(A3) 



g(t) = [i- 1 2 } 



d-3 1 



(A4) 



(t + iz) d - a+1 ' 
We want to show W a ^ when z — > 0. 
(i) d = 2m + 2 (m > 1) 

In this case g(t) has a branch cut on the real axis from —1 to 1. We perform the integration 
along the contour shown in Fig [10] (a), and obtain 



1 



dtg(t) = 7riRes t =- iz g(t) = ^ jjZT^Ai [jE^^ ~ 1 



d 



d—a 



d-3 



(A5) 



t=—iz 
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The derivative in (1A5j) can be calculated by the derivative of a composite function, 

|(d-a)] 



di<*- L f \ 2 )- r!( C /-a-2r)! (2tj ( lj 

(^)(^ - 1) " • • - (d - « - r - 1))(1 - t 2 )¥-(— ) 

(A6) 

where [|(d — a)] is the Gauss' symbol which is the greatest integer that is less than or equal 
to |(d — a). 

Then, when z — > 0, we obtain 

/ 1 , /n 1 -n d ~ a ,d — 3. ,d — 3 ,d — 3 ,d — a 

^(t) = --^(-l) ^ (_)(_-!)...(_ -(—-I)) 

d^a (A7) 

= -™(^y(4) 2 (^-3)(d-l) ■•■(«-!)• 

Then, from (jA"2l) . (TA~3]) and fTATD . for a = 21 < d(l e Z) is nonzero and W a has the form 
of (|22|) when z — > 0. (When a = d, we obtain dtg(t) = —ni from (1A5j) . (Note that g(t) 
has the branch cut, then [1 + z 2 ]( d ~ 3 )/ 2 — >■ (—1) when z — ► 0.)) 
(ii) d = 2m + 1 (m > 1) 

We perform the integration along the contour shown in Fig [10] (b), and obtain 

/ dtg(t) = -27riRes t= _ fa (/(*) - / <%(*). (A8) 
7-i Jc R 

For a = 2Z (Z £ Z), d — a is odd, so we obtain lim z _j,o Resj = _j 2i g(t) = from (1A6I) . Then, for 

a = 2Z and z->0we obtain 

/ dtg(t) = - [ dtg(t) = i(-l) d - a ^ d9e- l{d - a)e [l - e 2W }^ 
J-i Jc R Jo 

= i(-l) d - a (-2i) m - 1 [ d6e- i{2m+1 - a)e e l{m ~ 1)9 (sme) m - 1 



(A9) 



where 



-l) m 2 m - l i m {I sc [m - 1, m + 2 - a] - il ss [m - 1, m + 2 - a]) 

2 m - 1 i m+1 I ss [m - 1, m + 2 - a] ^ for odd m, 
2 m ~ 1 z m L fC [m — l,m + 2 — q]^0 for even m, 



/*7T /*7T 

I sc [m,n}= d6 (sin 9) m cos(n9) , I ss [m,n}= d9 (sin 9) m sin(n9). (A10) 
Jo Jo 



26 



Then, from (1A2j) (1A7j) and (1A9[) , for a = 21(1 G Z) W Q is nonzero and W a has the form of 
fl22|) when z 0. 

From (i) and (ii) we showed ( I2"2"j) for d > 3. Next we consider d = 2. 
(2) d = 2 

In this case we can perform the angular integral first, 



u •„(./'.//) — / / <iok { a e ikr 



(27T) 



^0 



i /-oo 1 /-oo 



(All) 



; 2 - a e- lk J (kr) = — / dxx 2 ~ a e~ zx J (x) 

'o 2vrH- a 7 

where Jo is the Bessel function of zeroth order. We perform the integral for a = 2 and 

a = 0. 

(i) a = 2 

In this case we have 

dxe- zx J (x) = . (A12) 



o 



Then, when z — > we obtain 



^ Q=2 (x, y) = W-^x, y) = -U (A13) 



(i) a = 

In this case we have 



°° dxxV 2 * J (x) = Q, 2, 1; -lj -> -1 (z 0) (A14) 

where F is the Gaussian hypergeometric function. Then, when z — >• we obtain 

W Q=0 ( Xj y) = W(x, y) = 2=^. (A15) 
Finally we have showed (1221) for d > 2. 



Appendix B: A formula for a finite series 

In this appendix we obtain a formula for a finite series by calculating the following integral. 

A= / dM; 2c / d^(sin^) 2 V fercose - £fe c>6>0 6,cGZ e,r>0 e,rGM. (Bl) 
Jo Jo 

This integral is a generalization of W Q in (lAlj) . The parameter e and r are auxiliary and they 
do not appear in the last formula. We obtain the finite series when we perform the 9 integral 
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before performing the k integral. On the other hand we obtain the simple expression when 
we perform the k integral before performing the 9 integral. Then we obtain the formula for 
the finite series. 

(i) We perform the 9 integral before performing the k integral. 
We perform the 9 integral, 

i d 2 r 1 d 2 

\2bikrcos9 l-y , L 11 \b jnikrcosB , L u \b 



d9(sm9Y b e tkrcos " = (1 + — — ) b \ d9e lkTCOSd = tt(1 + — — ) b J (kr) (B2) 

k z dr z J k dr z 

where Jo is the Bessel function of zeroth order. We substitute flB2j) into fIBlj) and perform 

the k integral. Then we obtain 

l=Q \ / JU l = Q \ / \ 

(B3) 

where fi = 2c — 21 + 1 and F is the Gaussian hypergeometric function. We have used the 
condition c > b > in the second equality in (1B3I) . When e — > 0, we obtain 

m F (a mi i- = 111,2 r ^ 



lim-^-F S^— — ,1; — - = ^ rz . (B4) 

e^o c m V2 2 ' ' e 2 / r(l - § )r(±±^) 



From ( 1B3|) and ( 1B4|) we obtain 

2« 



limA = 7r 3 / 2 y 6 Q r(/i) 1+ ( -) l = ^( 2c )!V k C, ^p— . (B5) 

(ii) We perform the k integral before performing the 9 integral. 
We change the variable as t = cos 9 and perform the k integral, 

"1 roo rl / i J \ 2c 



A = J' dt[l - t 2 }"- 1 / 2 dkk 2c e tkrt -* k = j 1 dt[l - t 2 }"- 1 / 2 ° J™ dke 



ikrt—ek 



j2c+l(^l f 1 dtll _ £ 2l6-l/2. 1 



r 2c+l L J (t + ^)2c+l 

(B6) 

where z = e/r. We perform the integration along the contour shown in Fig [TD] (a) in the 
same way as Eqs.( 1A5l) -( 1A7j) . and obtain 



Km A = ^ +1 ^t(--)^(-1) c (6 - \){b - |) . . . (b - c + f) 

= -^T^(& -|)(&- |) .-•(&- c + |) (B7) 

(2rV (-lY~ b 
= \i + i d I (26-l)!!(2c-26-l)!!. 
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From (1B5j) = (1B7I) we obtain the formula for the finite series. We simplify (1B5j) = (1B7j) 
and obtain the following formula, 



We can also rewrite (1B8I) as follows; 

B^-ff^W-^Sfef c^O^ (B 9) 
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' A-CB- X D CB- 1 
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